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We study superfluidity in a model of massless fermions coupled to a massive scalar field 
through a Yukawa interaction. Gap equations for a condensate with total spin J = are 
solved in the mean-field approximation. For the Yukawa interaction, the gaps for right- 
and left-handed fermions are equal in magnitude and opposite in sign, so that condensation 
occurs in the J p = + channel. At finite scalar mass, there are two different gaps for fermions 
of a given chirality, corresponding to condensation of particle pairs or of antiparticle pairs. 
These gaps become degenerate in the limit of infinite scalar mass. 

I. INTRODUCTION 

In fermionic matter which is sufficiently cold and dense, any attractive interaction at the Fermi surface 
leads to the formation of Cooper pairs Qj. The Cooper pairs form a Bose condensate, so that exciting a 
quasiparticle costs an amount of energy > 2\4>\, where |(^>| is the gap energy. This gap produces superfluidity, 
or, if the fermions are coupled to a gauge field, superconductivity. 

In quantum chromodynamics (QCD), one-gluon exchange between two quarks is attractive in the color- 
antitriplet channel. One therefore expects that at sufficiently large quark chemical potential /i and sufficiently 
small temperature T quarks condense into Cooper pairs which are color antitriplets. This condensate breaks 
the SU(3) color symmetry of the ground state, and gives rise to color superconductivity. For QCD, this 
phenomenon was first investigated by Barrois ||, by Bailin and Love j^], and others [Q. In seminal work, 
Baffin and Love estimated the gap energy to be on the order of \cf>\ ~ 10 _3 /i. Since the critical temperature 
for the onset of superconductivity, T c , is of the order of the gap energy, \<f>\, & color-superconducting phase 
of quark matter could form in the interior of neutron stars. Color superconductivity was also studied in an 
ST/ (2) gauge theory f§. 

Interest in this subject was renewed through works by Alford, Rajagopal, and Wilczek || and by Rapp, 
Schafer, Shuryak, and Velkovsky j?]]. In these papers, the attractive interaction was modeled with the 
instanton vertex. In the framework of a simple model of the type studied by Nambu and Jona-Lasinio (NJL) 
these authors found gap energies of the order of \(f>\ ~ /i. Since T c ~ \cf>\, for gap energies of this order of 
magnitude color superconductivity is no longer only relevant for astrophysical scenarios, it could also occur 
in relativistic heavy-ion collisions, whenever the bombarding energy is large enough to form relatively cold, 
baryon-rich quark-gluon matter. The maximum amount of baryon stopping was found in heavy-ion collisions 
at AGS energies, i?Lab ~ 10 AGeV; therefore, the most favorable conditions to form a color-superconducting 
phase of quark-gluon matter would occur in nuclear collisions somewhere in the energy range from i?Lab ~ 1 
to 40 AGeV, i.e., from GSI-SIS to CERN-SPS energies. 

Most studies of color superconductivity assumed that only the two lightest quark flavors, u and d, form 
Cooper pairs. This is based on the argument that, due to the larger strange quark mass, m s ~ 100 MeV, 
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the Fermi surfaces of non-strange and strange quarks do not match || , consequently there are no us or ds 
Cooper pairs. However, for a system where fi s = fi u = j-id = M, this mismatch is, to leading order in /z, only 
of order (kF, q — kp, s )l M ~ {m 2 s ~ m q)/^f j2 , and therefore vanishes at sufficiently large chemical potential. 
Based on this observation, Alford, Rajagopal, and Wilczek suggested a new form of color superconductivity, 
where all three flavors condense to lock the breaking of color and flavor symmetry ||. Comparison with 
earlier investigations jl(J of the patterns of symmetry breaking in chiral models revealed jllj] that the color- 
flavor locked state is energetically favored. This was independently confirmed by Schafer and Wilczek in a 
numerical study of the effective potential Jl^| . 

Other studies which investigated the interplay of color superconductivity and chiral symmetry breaking at 
finite temperature include those of Berges and Rajagopal [LJ and Langfeld and Rho H|. Color supercon- 
ductivity was also investigated via renormalization group techniques by Evans, Hsu, and Schwetz |l5| ] and 
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All of these studies were based exclusively on NJL-type models. The sole exception is the work of Son 
|l7j who applied renormalization group techniques to study the scale of the energy gap in QCD. He pointed 
out that non- instantaneous one-gluon exchange modifies the weak-coupling expression for the gap, \<p\ ~ 
H exp(— c/g 2 ) to \<j)\ ~ fi exp(— c/g); this was also noted in (TT]. This increases the likelihood to find color 
superconductivity at nonzero temperature. 

A quark-quark condensate A l j g aj3 is a N c x N c matrix in color space (i, j = 1, . . . , N c ), a Nf x Nf matrix in 
flavor space (/,<?= 1, ■ • ■ , Nf), and a 4 x 4 matrix in Dirac space (a, /3 = 1, ... ,4). Quark pair condensation 
can in principle occur in channels with arbitrary total spin J. Most previous studies indicate, however, that 
condensation in the channel with total spin J = is favored for two or more flavors. 

In this work, we shall not discuss the color or flavor structure of the gap matrix fig] , but rather focus 
on the Dirac structure. To do so, it is sufficient to study a simpler model, where fermions interact via the 
exchange of scalar bosons, 

C = 4>{i 1 -d-m)^-g^^+^(d^d»(j>-M 2 s <i ) 2 ) . (1) 

(We do not add a quartic self-interaction for the scalar field; such an interaction will not qualitatively change 
our results.) Scalar one-boson exchange is attractive, therefore, we expect the formation of Cooper pairs at 
finite density and sufficiently low temperature in the J = channel. Since the fermions in the Lagrangian 
([j]) are not charged, this model can only exhibit superfluidity, not superconductivity. Nevertheless, the Dirac 
structure of the superfluid condensate is closely analogous to that of the color-superconducting condensate 
in QCD. 

Our principal result is the following. Massless fermions have four types of J = condensates, corresponding 
to the pairing of fermions with the same helicity and chirality For fermions interacting via Yukawa 

interactions as in (|l|), the gap for right-handed, positive helicity quasiparticles is equal in magnitude (and 
opposite in sign) to that for left-handed, negative helicity quasiparticles. A novel feature of the relativistic 
treatment is the appearance of gaps for quasi-antiparticles. With Yukawa interactions, we find that the gap 
for right-handed, negative helicity quasi-antiparticles is equal in magnitude to that for left-handed, positive 
helicity quasi-antiparticles. These relations between the right- and left-handed gaps imply that there is 
condensation only in the J p = + channel. For M s < oo, the quasiparticle and quasi-antiparticle gaps are 
not equal. 

In the limit M s — > oo, the model ([j]) reduces to an NJL-type model, and the quasiparticle and quasi- 
antiparticle gaps become degenerate. Thus NJL-type models are unrealistic, in that they force the equality 
of the quasiparticle and quasi-antiparticle gaps. 

The number of condensates is not merely a technical matter, but can even affect the order of a su- 
perconducting phase transition. Consider the transition where up and down quarks condense to a color- 
superconducting phase. In we demonstrated that if only one condensate appears, such as that for 
particles, the transition is in the universality class of a single C/(3) vector and can be of second order. Based 
on the analysis of this paper, it seems unavoidable that condensates for both particles and antiparticles 
appear together. The universality class becomes that of two U(3) vectors; while the fixed-point structure 
of this model is unknown, it might be driven first order by the Coleman- Weinberg phenomenon. If true, it 
implies that for two degenerate flavors, the color-superconducting transition is of first order, independent of 
the chiral transition |L3| . 

In [ pH we also showed that the transition where up and down quarks condense with strange quarks is 
driven first order by the Coleman- Weinberg mechanism. Schafer and Wilczek argued that the degrees of 
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freedom in quark and hadron matter match at nonzero /i [fl2[ . Their observation implies that this line of first 
order transitions, which starts at zero temperature, could terminate at some nonzero temperature, similar 
to the liquid-gas transition in water or nuclear matter. 

The outline of the paper is as follows. In Section II we compute the full fcrmion propagator and discuss the 
excitation spectrum in a superfluid. Then, the gap equations are derived in the mean-field approximation 
for the model (Q). We solve them in Section III. A logarithmic ultraviolet singularity arises, which upon 
renormalization introduces a renormalization scale A. We study solutions of the gap equations as a function 
of the coupling g and the renormalization scale A. We also analyze the temperature dependence of the gap 
energy and discuss the critical temperature T c , at which superfluidity vanishes. Section IV concludes with 
a summary of our results. Appendix A contains a derivation of the grand partition function and the gap 
equation in the mean-field approximation, both for an antifermion-fcrmion condensate and a fcrmion-fcrmion 
condensate. In Appendix B we analyze the Dirac structure of the gap matrix in detail and make contact 
to earlier work In Appendix C we discuss the renormalization of the gap equations. Our units are 
% = kg = c= 1. The metric tensor is g^ = (+, — , — , — ). 



II. THE GAP EQUATIONS IN THE ULTRARELATIVISTIC LIMIT 

A standard way to study superfluidity is to compute the grand partition function Z(T, V, fi) associated with 
lh in the mean -field approximation for a fermion-fermion condensate. This is done explicitly in Appendix 
\, cf. cq. ( A15 ). One arrives at an effective action of the form (cf. eq. ( A16 ), see also eq. (1.1) of g]) 

=[ ${ x ) (*7- 9 + ^7o - m)ip(x) + - / [ipc(x)A(x,y)ip(y) +h.c] , (2) 

J X J x.y 

where ipc(%) is the charge-conjugate spinor, defined by 

1>c{x) = C$ T {x) , i > c{x) = ij?{x)C , 4,(x)=C4%(x) , $(x) = 1%(x)C , (3) 

and C — i7 2 7o is the charge-conjugation matrix in Dirac representation, C = — C _1 = — C T = — C>. We 
also abbreviated J = d(it) J v d 3 x. A(x, y) is the gap matrix. If the system is translationally invariant, 
A(x,y) = A(x — y), the Fourier transform 



A(fc) = / e lk ' x A(x) (4) 

J X 

obeys the gap equation (cf. eq. ( A35| ), see also eq. (1-34) of ||) 

A ( fc ) = 9 2 1 E D( - k - M *) G o (l) A (9) G + (l) ■ (5) 

Here, fc = -t(2n + 1)ttT, J2k = En V J d 3 k/(27r) 3 , and 

[ G o] ( k ) = 7 ' fc ± M7o - m (6) 
is the free inverse fermion propagator, while 

^{[Gf]- 1 -^}" 1 (7) 

is the fermion propagator dressed by the interaction with the fermion-fermion condensate, E = 
A ± , A + = A, A - = 7oAt7o- Gg~ corresponds to propagation of free particles, G$ to that of 
charge-conjugate particles. G + propagates quasiparticles and charge-conjugate quasiparticles, respec- 
tively D(p, M s ) = l/(Mg — p 2 ) is the propagator of the scalar boson. 

In Appendix [b| we show that massive fermions have eight possible gaps in the J = channel. In the 
ultrarelativistic limit, these reduce to four. They correspond to the pairing of right-handed fields with 
positive helicity with themselves, etc. Consequently, the gap matrix (0) has the form 
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A(fc) = 0+ + (k) V+ + (k) + 0+ (*) 7>+_ (k) + ^_ (k) P~_ (k) + (*) P." (k) 



(8) 



Here, 

P+ + (k)=7> r 7>+(k) , P+_(k)=W_(k) , P"_(k)=P r P_(k) , p- + (k)=P^P + (k) , 
are projectors onto states with given chirality and helicity, as denoted by the two subscripts, where 



V r ,i = 



1±75 



p±(k) = 



1 ± 757o7 ' k 



(9) 



(10) 



The additional superscript in (|^) refers to the sign of the energy for a non-interacting fermion field. For 
massless fermions, this superscript is in principle superfluous: right- (left-) handed fermions with positive 
(negative) helicity must have positive energy, while those with negative (positive) helicity must have negative 
energy. Nevertheless, we find it physically illuminating to retain it. 

The quantities 4>fe± are the individual gap functions; for example, (f>^ + corresponds to the condensation 
of two right-handed particles with positive helicity, while <p~_ corresponds to the condensation of two right- 
handed antiparticles with negative helicity. This means that in the J = channel, fermions of different 
helicity or chirality do not condense. This is due to the projectors in eq. (H): the effective action (0) reduces 



to a sum of four terms, one for each given chirality and helicity, sec cq. (B34). 

The fact that only fermions with the same helicity condense in the J = channel is physically obvious in 
the center-of-momentum frame of the pair, and holds for either massive or massless fermions. In that frame, 
the requirement of equal helicity means that the projections of the spin along the direction of motion are 
antiparallel, as in the usual non-relativistic treatment of superconductivity. Taking the fermions to move in 
the z direction, in the center-of-momentum frame L z — 0; thus if J = 0, also S z — 0. The spin wave function 
with S z = can be either the antisymmetric singlet, S = 0, or the symmetric triplet, 5 = 1. To obtain total 
spin J = 0, for the former, S = combines with L = 0, while for the latter, S = 1 combines with L = 1. 
The signal for L = 1 is the appearance of k in eqs. (|J), ([To|). 

In the following, we first discuss the structure of the full fermion propagator (0), because the poles of the 
full propagator determine the excitation spectrum in a superfiuid. We then return to the gap equation (||). 



A. The full fermion propagator 

In the ultrarelativistic limit, the full propagator (^) assumes the form 

G + (q) = 7 ■ q + ^7o - 7o [A(q)] f 7o (7 ' 9 - A^o)" 1 A(q) 
From eq. (0) one computes 



7o At 7o = [0+ + ] t ^T f +[0+] t 



V~ + \d>-_] Vf 



(11) 



(12) 



Since the condensates <p are not matrices (i.e., unlike the QCD case, they do not carry other internal degrees 
of freedom like color and flavor), the hermitean conjugation is replaced by simple complex conjugation. With 
the identities 



and 



7 oAt 7o ( 7 . q _ M7q) A = (7 • g — M7o ) (|0+ | 2 V+ + + \4>j_\ 2 V+_ + |#T_ | 2 Vf_ + |^ + 1 2 Vj + ) (13) 



(7 • 9 ~ M7o) (7 • q + M7o) = [qo - (|q| - A*) 2 ] {V+ + + V+_) + [<? 2 - (|q| + /i) 2 ] {K- + K+) > ( 14 ) 



one derives 
G+(q) 

where 



^ + +(q) 



KM) 



*t)]' 



(0r-)] 9o - [e {<t> t+ j\ 



+ 



^+(q) 



- Nl2 



(7 -q- nio) , (15) 
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f 



= [(|q| TM ) 2 + |0| 2 ] 1/2 . (16) 



The full propagator ( pL 5| ) has eight poles in the complex qo plane: 

q^=±e+ (<(>++) , qi 3 ' i) =±e+(ct>t) , qi 5 ' 6) =±e-(4>-_) , q™ = ±<T(#- + ) . (17) 

These poles deserve further discussion. The kinetic energy of non-interacting, massless particles with 3- 
momentum q is |q|. Let w(|q|) denote the energy of the quasiparticle excitations as a function of the kinetic 
energy of the non-interacting particles. In a non-interacting, ultrarelativistic Fermi system quasiparticles 
(quasi-antiparticles) are the usual particles (antiparticles) and consequently the excitation spectrum is 

tL> (|q|)— =+|<l|— M particles, 

= — 111 + A* particle holes , , . 

= — | q| — antiparticles , ^ ' 

= + IqI + A 1 antiparticle holes . 

On the other hand, in a superfluid the pairing of particles and antiparticles changes the excitation spec- 
trum. As is well-known the excitations are now quasiparticles (quasi-antiparticles) which are in essence 
a linear superposition of particles (antiparticles) and particle holes (antiparticle holes). This linear superpo- 
sition is the well-known Bogoliubov transformation jlj. The quasiparticle excitations are sometimes called 
"Bogoliubons" . On account of ( |l7| ) , in a superfluid the spectrum of right-handed excitations is 



w(|q|) — fi = — e + ((f>r + ) = — y (|q| — ^) 2 + |</V+I 2 right-handed quasiparticles , 

= + e + (0+ + ) = + J (|q| — /i) 2 + |</>.^ + | 2 right-handed quasiparticle holes , 

= — e~((f>~_) = — J (|q| + fj,) 2 + \(/>r-\ 2 right-handed quasi-antiparticles , 



(19) 



= + e~(<fi~_) = + y (|q| + A*) 2 + |0^-| 2 right-handed quasi-antiparticle holes 

while for left-handed excitations one has 

w(|q|) — A* = — e+ (0/-) = ~ \J (|qj — m) 2 + I0^L| 2 left-handed quasiparticles , 

C +(A& \ — -j- . /(|q| _ ^2 _|_ 1 | 2 left-handed quasiparticle holes , 



= — e (4> l+ ) = - y (|q| + A 4 ) 2 + l^+l 2 left-handed quasi-antiparticles , 

= +e~(<fij + ) = + \ (|q| + m) 2 + l<fc-l 2 left-handed quasi-antiparticle holes 



(20) 



These branches are displayed in Fig. [I], for simplicity with a common gap 4> — 0.5 /i. For free fermions, 
the particle and hole branches cross at the Fermi surface. The generation of a gap produces a phenomenon 
analogous to level repulsion in quantum mechanics, as the quasiparticle and quasiparticle-hole branches 
become disconnected. Relative to the non-interacting case, the meaning of quasiparticle and quasiparticle- 
hole branch interchanges above the Fermi surface. 

The main difference from the non-relativistic case is the appearance of the (quasi-) antiparticle branch 
and the (quasi-)antiparticle-hole branch. This means that condensation is not only restricted to particles 
close to the Fermi surface, but that also antiparticles in the Dirac sea condense. This can be confirmed by 
considering the fermion-fermion scattering amplitude: a BCS-type singularity occurs not only for particle- 
particle scattering, but also for antiparticle-antiparticle scattering, giving rise to the antiparticle condensates 
<t>r- £+■ Notice, however, that it is always easier to excite a quasiparticle than a quasi-antiparticle. It costs 

only an energy 2 </> to excite a quasiparticle at the Fermi surface, but 2\J /i 2 + (j) 2 to excite a quasi-antiparticle 
with zero momentum. 
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FIG. 1. The dispersion relation for (quasi)particle (solid), (quasi-)antiparticle (dotted), (quasi)particle-hole 
(dashed), and (quasi-)antiparticle-hole states (dot-dashed), (a) Non-interacting, massless fermions, (b) superfluid, 
massless fermions. 

One can also compute the occupation numbers for quasiparticles and quasi-antiparticles. To this end, it is 
advantageous to Fourier-transform the full propagator to the so-called "mixed" (or "Saclay" ) representation 



G+(r,q)=T^e^G+K,q) , 



(21) 



where iu n ~ i(2n + 1)t:T = —qo- For instance, the free propagator for massless fermions reads in this 
representation |19[ 



G +(r,q) 



9{r) 



N F 

-N F 



T 

|q| 



T 



cxp[-(|q| -v)t] A+(q) 7o 
exp[(|q| +v)t] A"(q) 7o , 



(22) 



where Np(x) = (e x + l)^ 1 is the u sual Fermi-Dirac distribution and A ± (q) is the projector onto positive 
(negative) energy states (eq. ( |B22 ) in the massless limit, i.e., for /3 q = 1 and a q = 0). The occupation 
numbers are now directly exhibited as the prefactors of the exponential (imaginary) time evolution factors. 
For r > 0, the first term corresponds to the propagation of particle holes, with the thermal occupation 
number 1 — A^[(|q| — /i)/T], while the second term corresponds to the propagation of antiparticles, with the 
thermal occupation number A^[(|q| + fi)/T]. For r < 0, the first term corresponds to the propagation of 
particles, with the thermal occupation number ./Vp[(|q| — fi)/T]. The second corresponds to the propagation 
of antiparticle holes, with occupation number 1 — Np [(|q| + /x) /T] . 

The full fermion propagator is computed from eq. (|2l|) with (|15|). The result is 



G+(r,q) = - £ 

hs=r+,l- 



E 

hs=r-,l+ 



(r) - N F 

-r) - N F 
(t) - Np 



T 



[l-n+(0+)]exp [-e^ 



J hs> 



T 



<(4>&) ex P [ e + (^)r]^+(q) 7 o 



T 



[1 



,)] ex P [" 



G 



where 



(-t) - N f 



,J hs) 



T 



n^hs) ex P [ e i < t ) hs) T ] ) ^ hfl (q)7o 



e ± WT(|q|Ti") 



2e±(0) 

In deriving eq. (^3|), in an intermediate step one encounters terms of the form 

(e ± -M)7o-7-q , -(e ± +M)7o-7-q 
2^ and ^ ■ 

By adding and subtracting a term |q| 70 in the numerator, one can transform them into 



(23) 



(24) 



7o 



l-n±)TA^(q) 7o ^ 



and 



- 7 o< T AT(q) 7 o^ 



Eq. ( p3|) then follows by making use of eq. ( B31 ) and the orthogonality of the energy projectors, A + A~ = 0. 

There is a one-to-one correspondence between the eight poles of the full propagator ([l7]) and the eight 
terms on the right-hand side of eq. (23). The first line in eq. ( p3|) arises from the pole for right (left-- 
handed quasiparticle holes, the second line from that for right(left)-handed quasiparticles, while the third 
line originates from the pole for right (left)-handed quasi-antiparticle holes and the fourth line from that for 
right(left)-handed quasi-antiparticles. As in the non-interacting case (p2[), the occupation numbers can now 
be directly read off as the prefactors of the exponential time evolution factors. 

It is interesting to study two limiting cases of (|2^)- The first is the zero-temperature limit, where the Fermi- 
Dirac distributions vanish on account of the fact that their argument is always positive. Then, the functions 



defined in eq. 



are identified with the occupation numbers for quasiparticles and quasi-antiparticles 



in a superfluid at T = 0. Correspondingly, 1 — are the occupation numbers of quasiparticle holes 
and quasi-antiparticle holes, respectively. These are shown in Fig. |2|. The quasiparticle and quasiparticle- 
hole occupation numbers exhibit the smearing around the Fermi surface characteristic for a superfluid or a 
superconductor jjj . The smearing is a consequence of the fact that the "Bogoliubons" are superpositions of 
particle and hole states. Note also that there are always some quasi-antiparticle-hole excitations present. 
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FIG. 2. The occupation number for quasiparticles (solid), quasi-antiparticles n q (dotted), quasiparticle holes 
1 — rit (dashed), and quasi-antiparticle holes 1 — n~ (dot-dashed) in a superfluid at T = 0. 
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The other interesting limit is the non-interacting case where one is supposed to recover eq. (|2 
all gaps to vanish, (fy^ s — 0, one has e ± — > ||q| =p fx\, such that 



n+-ffl(M-|q[ 



1 (0-0) 



Taking 
(25) 



The contribution from quasi-antiparticle holes vanishes completely, and the quasi- antiparticle term becomes 
the antiparticle term, the second line of eq. (22). However, in order to obtain the particle-hole contribution 
(the first line of eq. (|2^)), one has to combine the contribution from quasiparticles and quasiparticle holes. 
This can be understood by inverting the Bogoliubov transformation which combines particles and particle 
holes to quasiparticles (and quasiparticle holes). 



B. The gap equations 

After discussing the full fermion propagator and the excitation spectrum in a superfluid, we now return 
to the discussion of the gap equation (||) in the ultrarelativistic limit. With the expression ([l5|) for the full 
propagator, it assumes the form 



A(k) 



,T 
V 



Y,D{k-q,M s ) 



<#+(<?) 



«8 



[e+Wv + +)] 
4>r- (q) 



4 



2^+(q) 



#-(«) 



(q) 



j2?V-(q) 
T2^ + +(q) 



95 - L e ~ 



(26) 



Forming suitable projections of eq. (||), one derives gap equations for the individual gap functions: 
9 2 T 



<p;_(k) 
07 + (fc) 



2 V 

2 V 
9 2 T 



tl 
2 V 



q,M s ) 
J2 D(k-q,M s ) 
q,M s ) 
q,M s ) 





i k 


q 








1 + k 


q 


q 2 








ql 




V + )] 2 




i -k 


q 






h — 


1 + k 


q 


q% 


- 


++)]' 






- [e-(< 


K-)] 2 




1 + k 






(<zH 




i - k 


q 


q 2 


- 








- 






1 + k 


q 




.(?)- 


1- — 


1-k 


q 


4 




V + + )] 2 








*r-)] 2 



0r- (<?) 

4>r- (q) 



(27a) 
(27b) 
(27c) 
(27d) 



These gap equations (^) do not agree with those derived by Bailin and Love || . These authors include 
only half of the fermion quasiparticle modes: they take the quasiparticle and quasiparticle- hole branches, 
with gaps <j)^ + and <j)^_ , but neglect the quasi-antiparticle and quasi-antiparticle-hole branches, with gaps 
<p~_ and ■ They also restrict themselves to condensation in the + channel, which forces (f>^ + — —<j>j_ 
(see eq. ( B35| )). Doing so, they find a single equation for a linear combination of the gaps. We do not find a 
s ing le equation. Instead, due to the form of the full fermion propagator (|l5|), each of the four gap equations 
(|27j) contains a sum of two terms: each gap on the right-hand side, such as (f'r+i IS always multiplied by an 

energy denominator with only that gap, 1/ |(/q — [e + (0^ + )] 2 |. 

In the case of a point-like four-fermion interaction as in NJL-type models, D(k — q, M s ) — l/M^. Then, 
the terms proportional to k • q vanish by symmetry, and the right-hand sides of ( [27] ) no longer depend on 
either fco or k: the gaps are simply constants. Defining 



n<t>) 

the gap equations reduce to 



T 



2M ? V 



E 



ql 



(28) 



S 



0+_ = </>;_ = 4>t F(<Pt) , 4>t = 4>J+ = <t>t+ H<Pt+) ■ (29) 

The gaps for right- and left-handed quasi-antiparticles equal those for the corresponding quasiparticles. This 
reduces the number of independent gaps to two. 

The solutions of eqs. ( p9| ) satisfy 1 = F(4>)F(4>F((f))), where 4> is either 0+ + or 4>f_. One possible solution, 
4>o, satisfies 1 = |F(0o)h as F(4>) is single-valued as a function of <p for <f> > 0, and F(6) — F (—(/>), the only 
possible solutions are ±4>o, therefore 4>^ + = ±<l>f_- Since F < 0, the solution of eqs. ( |29| ) obeys (f>f_ = — 
From eq. (B35) we then conclude that condensation occurs only in the + channel. 



In conclusion, for scalar NJL-type interactions there is only one independent gap function. Physically, 
this happens because for a point-like interaction there are no states with nonzero angular momentum, and 
so terms in the gap proportional to k — which signal L = 1 — must vanish. In the next section we show 
that scalar boson exchange over a finite range, M s < oo, lifts this degeneracy and produces two independent 
gap functions. 



III. SOLVING THE GAP EQUATIONS 



The set of equations ( P7|) determines the gap functions for massless fermions in the mean-field approxi- 
mation. In principle, the gaps are functions of the 4-momentum A:^, so that eqs. ( p7| ) are actually integral 
equations. In the following, we assume that the momentum dependence of the gap functions is negligible. 
We comment on the validity of this approximation below. Remember that in NJL-type models the gap 
functions are automatically independent of fc M . 

Under the assumption that the gap functions are simply constants, one ends up with the following set of 
four algebraic equations: 



t- [^o 

r+ [Fq 



J r+ 



i+J-^ + (0t)]+^+ 

b+ + ) -Ft {&+)]+ fc- 



-)+-7T( 
-)--7T( 



where 



T 
V 



(k-q) 



*7 + )] 
4>r-)] 

ft-)] 



(30a) 
(30b) 
(30c) 
(30d) 

(31) 



Note that = for NJL-type models. 

The functions are further evaluated replacing the Matsubara sum over q$ = —i(2n + l)nT by a 
contour integral and applying the residue theorem. Poles in the complex q plane arise from the fermion as 
well as the boson propagator. Keeping only the former, we obtain 



J-Q — 



9 



Ft 



327T 2 fc 

9 2 

327T 2 fc 



dq In 



Ml + {k + q) 2 



M 2 + (k- qf 



-j: tanh 


\ 6± ] 


e ± 


2T 



dq 



M 2 



2qk 



In 



M 2 S + (k + qf 
M 2 +(k- qf 



—t; tanh 

€ ± 



2T 



(32a) 
(32b) 



where we have performed the angular integration and denoted k = |k|, q = |q|. Since JF 1 are real, the gap 
functions 4> can be chosen to be real. 

To obtain this result we assume that the exchanged boson has zero energy, qo = fco- This approximation 
can be justified as follows. The dominant term in Jqi comes from particles close to the Fermi surface, q ~ /i. 
Assuming that \<f>\ -c fJt, this produces a logarithmic dependence on \<p\: 



dq 



dq 



In 



A 



uv 



\<P\ 



(33) 



(We discuss the ultraviolet cut-off Ayv A* below.) The assumption that the gap \<p\ is small relative to 
the Fermi energy is justified in weak coupling, g <C 1. To obtain the logarithm d33|), the fermions have to 
stay close to the Fermi surface. This can only be achieved if the exchanged boson has zero energy, go = 
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Poles of the boson propagator, q — k ± \fM~i + (k — q) 2 , which we neglected in deriving eqs. (p^) , 
represent processes in which the fermions in the gap equation are far from the Fermi surface; this produces 
terms of order 1, but not of order ln(Auv/| </>!)■ (Further, keeping only zero-energy bosons also eliminates a 
possible dependence of on k .) 

Analogous to eq. (|33|), the functions Tq X behave as 



dq 

e~ 



dq 



In 



A 



uv 



(34) 



in the limit Auv ^S> ^ ^S> \tj>\- Therefore, !F 01 do not include terms ~ ln(Auv/|<^|) and, in weak coupling, 
can be neglected relative to Physically, this is because e~ represents the excitation spectrum of quasi- 

antiparticles which are always far from the Fermi surface, cf. Fig. 

While our approximations are controlled only in weak coupling, we nevertheless find it illustrative to 
consider the qualitative nature of the solutions in strong coupling. For g ~ 1, the gap can be of order \x. 
Similarly, the functions Tq 1 are of comparable magnitude to . We therefore retain them in the following 
analysis. 

It is surprising that an ultraviolet cut-off Auv appears in the gap equations. There is no ultraviolet 
divergence for the gap equations of ordinary superconductors JD], as a cut-off is provided by the Debye 
frequency, lod <C fi, so that the integrals receive contributions only from a narrow interval around the Fermi 
surface, /i — lud < q < /i + lujj. On the other hand, in He-3 the cut-off is provided by the chemical potential, 

The appearance of Auv is an artefact of our approximation of taking a constant gap. This approximation is 
manifestly inconsistent when k 3> /i. The gap functions fall off at large k, removing the apparent ultraviolet 
divergence. Consequently, the true gap function is proportional to /i, not Auv- That there is no ultraviolet 
divergence in the true gap function can also be seen by considering the particle-particle and antiparticle- 
antiparticle scattering amplitudes. The box diagrams whose infrared singularities generate a non-zero gap 
jjj are manifestly ultraviolet finite, even in the relativistic regime. In order to simplify the solution of the 
gap equations, however, we take the gap functions to be constant and thus ignore these complications. 

With this approximation, the integral is logarithmically divergent in the ultraviolet, while is finite. 
Instead of using a fixed ultraviolet cut-off Auvj i n Appendix ^| we show how J 7 ^ can be rendered finite by 
renormalization. The result is 



J n 



9 



32n 2 k 



dq 



In 



M 2 +(k + q) 



2 1 



Mi + (k- qf 





f C± l 


tanh 






2T 



In 



{k + qf 
{k-qf 



16tt 2 A 2 e 2 ' 



(35) 



where A is a renormalization scale. One can readily convince oneself that for A ^> fc, the renormalization 
scale can be identified with the ultraviolet cut-off Auv introduced above. In light of this, we take A ^> [i. 

In the following, we turn to the solution of the gap equations, first at T = and then at nonzero T . It is 
instructive to start with the weak-coupling limit, g <C 1. 



A. Weak-coupling limit 



The gap energy <p is expected to be exponentially small in weak coupling, (f> ~ //exp(— c/g 2 ) Jl|. The 
integrands of the functions Jqi are strongly peaked around q = fj,, since then e + = \<p\ <C /i. Consequently, 
the main contribution to the integrals comes from a (small) region n~5<q</i + S, where S = a/i, 
with some constant a which we do not determine. The functions are relatively suppressed by a factor 



*+/e- ~ |<£|/2/i< 1. For<5» \<j)\: 



•> n 



9 M 



In 



Ml + {k + n) 2 
M s 2 + (fc-/x) 2 



16tt 2 k 



g 2 ii f M 2 + k 2 + fi 2 







16tt 2 k 
T7 



2fj,k 



In 







1 26 

ln H' 

M 2 + (fc + / i) 2 



M 2 + (fc- M ) 2 



2(5 



(36a) 

(36b) 
(36c) 



Renormalization corrections are unimportant in this limit, since they only change 6. 
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The gap equations ( pOj) were derived under the approximation that the gaps are constants independent of 
k. From eq. (|3^), however, we see that the functions J-qi do depend strongly on k, and peak around k = \i. 
Thus, as was seen previously in eq. ([33]), in weak coupling pairing is dominated by fermions close to the 
Fermi surface. For g <C 1 it is therefore consistent to neglect the momentum dependence of both !Fqi and 
the gap functions and consider the above expressions (|3^) at k = fi. 

To one-loop order, a scalar mass is generated by its coupling to a fermion loop, 



M 2 = 9 -(^ 
2 U 2 



T- 



(37) 



If we had included scalar quartic interactions in the Lagrangian ([!]), they would contribute to M 2 a term 
~ AT 2 , where A is the quartic coupling. Thus, at zero temperature there is no change, while at nonzero 
temperature there is an innocuous increase of M 2 . As will be disussed below, this tends to decrease T c . 
In weak coupling, M s ~ <?/i <C /x, and we can further approximate 



<r 



•> l 



g 2 , V 2 2<5 

o In — - in — 

16tt 2 M? 



1*1 



We can now solve the gap equations 



16tt 2 
i. Using 
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In- 



u = - 
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8vr 2 




Pr+ 


9 2 
8tt 2 


r- = - 


Pt 


9 2 
8tt 2 


>e+ = - 
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9 2 
8tt 2 



In 



In 



!^) we obtain 
26 



26 

W\ 



\<pt 

26 



In 



In 



1^ 
4 M 



M| 
M 2 



In 



- 1 



In- 



2<5 



10? 



(38a) 
(38b) 

(39a) 
(39b) 
(39c) 
(39d) 



By the same arguments used at the end of section O in the case of a point-like four-fermion interaction, 
the first two equations yield = —</>£_, where ^ + is a solution of 1 = g 2 /(8n 2 ) ln(2£/|0i_|). The 
last two equations indicate that the gap functions for the quasi-antiparticles are larger than those for the 
quasiparticles by a factor ln[4// 2 /M 2 ] — 1. In conclusion, the solution of the gap equations to leading order 
in weak coupling reads: 





" 8tt 2 " 




" 8^ 2 " 


exp 




~ /i exp 






. 9 2 _ 




. 9 2 . 



i 4 ^ 2 i 



(40a) 
(40b) 



By eq. (B35), this confirms the result already obtained in the limit M s — *■ oo that J = pairing of fermions 
interacting via scalar boson exchange only occurs in the + channel. Note, however, that in contrast to the 



M s —> oo case, where 
4>~_ - ln(l/V^ + 



J r+ : 



in weak coupling there are two independent gaps instead of one, with 



B. Strong coupling 

In strong coupling, fermions of all energies |q| contribute to the functions so that they have to 

be computed numerically. Since these functions depend on k = |k|, the set of algebraic equations (|30| ) is 
consistent only if this k dependence is negligible. Before actually solving (pTJ), we therefore investigate the 
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functional dependence of !Fq X on k in detail. This dependence is shown in Fig. || for g = 10 and in Fig. [| 
for g — 1 for various values of <p. The renormalization scale is A = 10 /i in both cases. 
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FIG. 3. The functions as functions of k/[i for (p = fj, (solid), 0.1 /i (dotted), and 0.01 (dashed). The coupling 
constant is g — 10, the renormalization scale A = 10 n, the temperature is taken to be T = 0. 
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FIG. 4. Same as in Fig. but for <? = 1. 

As the coupling g or 4> become larger, the functions ^ depend less strongly on k. In weak coupling 
and for small cj>, remain approximately constant, although they are then small in magnitude, while 
^0~j become peaked around the Fermi surface. In strong coupling, is smaller than since then 
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M s ~ gfi becomes large and the theory approaches the NJL limit, where JF-f = 0. Note that changes in the 
renormalization scale shift the values of the functions J-^ by a constant amount; a larger A increases the 
absolute values of these functions. 

Lastly we find that while the functions change with 0, the functions T^ x are nearly independent 
of <fi. The latter is easily understood for small \<fi\ <C (J>, since only enters the integrals through er = 
\/(q + n) 2 + <j> 2 — 2/x for \cj)\ <C /i. It is interesting, however, that this behavior persists also for \<j>\ of order 
fi. 

In conclusion, neglecting a possible k dependence does not appear to be a bad approximation unless 



the coupling becomes substantially smaller than 1. The results of subsection III A are then sufficient to 



obtain the exponential dependence of the gaps on the coupling; to determine the prefactor S requires more 
sophisticated methods to solve the gap equations, for instance the approach of Ref. pc| ]. However, the gap 
energies become rather small for g < 1, cf. Fig. |[ We shall therefore numerically evaluate the gap equations 
( |30| ) only for g > 1, where the k dependence can be safely neglected. In what follows, we then always take 
the values of the functions J-^ l at momentum k = /i. 

Solving the gap equations ( |30| ) numerically, we find that the solutions satisfy the relations 
as is true in the limits of M s 



r±_ 



—<j>\_, cf)~_ = —0^, as is true in the limits of M s — > oo and weak coupling. According to eq. ( B35 ), 
this implies that for massless fermions interacting via scalar boson exchange, condensation is possible only 
in the + channel. 



In Fig. U (a) we show the dependence of the quantity q 

on g 2 /47r, and in Fig. [s] (b) the associated behavior of 
values of the renormalization scale A. 
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FIG. 5. (a) The dependence of (j> (in units of fi) on g 2 /4ir for A = 10 fx (solid) and 100 fj, (dotted), (b) The 
corresponding values of cj>~_ (upper set of curves) and 4>t+ (lower set of curves), in units of <f>. The dashed curve in 
(a) corresponds to the weak-coupling limit for <f>. 

The value of (j) decreases rapidly with the coupling, but values of <j) of the order of \i and larger are possible 
if the coupling is sufficiently large. The value of <p increases with increasing renormalization scale A. The 
overall behavior is in good agreement with the weak-coupling limit 



= 2(W2 + 21n 



4/i 2 



M? 



1/2 



exp 



8^ 
9 2 



(41) 



which is shown by the dashed line. Here, S ~ 25 \i was taken to fit the numerical results for A = 10 fi. For a 
respectively larger value, the weak-coupling limit can also approximately fit the results for A = 100 /i. 
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Turning to Fig. || (b) we observe that, for increasing coupling, the difference between <f>~_ and </> + + 
decreases. This is in accord with our expectations, since the mass increases ~ g, and consequently, the 
boson propagator approaches the form D(p, M s ) ~ 1/M| . We have already discussed at the end of section 
H that in this case <p^ + = (f>~_ and <j>^_ = <j)J + - An interesting observation is that the values for the quasi- 
antiparticle gap <pZ are larger than those for the quasiparticle gap </>+ + for all values of g, not just in weak 
coupling, cf. eq. (pOj) . 



C. Temperature dependence of the gap 



Finally, we discuss the temperature dependence of the gap functions. As in ordinary superconductors, 
the value of (j> decreases with temperature, and vanishes at T c , the critical temperature for the onset of 
superconductivity (or, in our case, superfluidity). T c can be estimated in weak coupling, where the integrals 
peak in a narrow region around /i. With cf> = one derives 



•> 1 



g 2 , V cs 
16^ ln MF ln i; 



167T 2 



In 







a^ 1 

Ml 
. 



hi 



(42a) 

(42b) 
(42c) 



where Q = 2e 7 /it, and 7 is the Euler-Mascheroni constant. Using <p^ + = —4>t_, the gap equations at T c yield 
the condition 



1 



9 2 , (S 
— 77 In — 

SlT 2 T c 



for the critical temperature, i.e., 



0+ (T = 0) -0.570+ (T = 0) 



(43) 



(44) 



This relation is identical to that found in BCS theory H. 

For arbitrary the value of T c has to be determined numerically. For the temperature-dependent boson 
mass considered here, the increase of the mass with T actually leads to smaller values for T c than expected 
from the BCS result (phf), see Fig. || This is in agreement with the general expectation that larger boson 
masses lead to smaller values for the gap 0, and consequently to smaller values of T c . We find that and 
4>~_ always vanish at the same temperature for any value of the coupling. 
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FIG. 6. The temperature dependence of 4>t+ an d 4> r - for g = 5, A = 10 /x. T c ~ 0.47 fi is smaller than the value 
0.57 (j)^ + (T = 0) ~ 1.25/1 expected in weak coupling. 

IV. CONCLUSIONS 

In this paper, we have investigated superfluidity in a system of massless fermions interacting via scalar 
boson exchange. The gap matrix contains in general four independent gap functions, corresponding to the 
condensation of fermions with the same helicity and chirality. We solved the gap equations in the mean-field 
approximation as functions of the coupling g and the renormalization scale A. For scalar boson exchange, 
condensation in the - channel does not occur, and the number of independent gap functions reduces to 
two, one for quasiparticles and one for quasi-antiparticles. The quasi-antiparticle gap is found to be larger 
than the quasiparticle gap, by a factor ~ ln(l/g 2 ) in weak coupling. This is in contrast to NJL-type models, 
where the point-like four-fermion interactions do not allow for pairing with L = S = 1 and force the equality 
of the quasiparticle and quasi-antiparticle gaps. We also analyzed the temperature dependence of the gap 
functions, and found that a temperature-dependent boson mass can significantly reduce T c . 

Qualitatively, our model shares the feature with QCD that the gap is exponentially small for g <C 1. 
As can be seen from Fig. || (a), in that regime small variations of the coupling lead to order-of- magnitude 
changes in the gap. Moreover, the gap energies are quite sensitive to the value of the ultraviolet cut-off 
Atjv ~ A. Taking a small value for the coupling, Bailin and Love found gap energies which are extremely 
small, ~ 10 -3 ^ ~ 1 MeV for fi ~ 1 GeV (3|. In more recent studies j6|0||l3|], much larger gaps, ~ 100 
MeV, were found. Such large gaps can be obtained from the treatment of Bailin and Love simply by taking 
larger coupling constants. Alternatively, the latter studies would get smaller values for the gap by reducing 
the coupling constant, or changing the form factors employed there to remove the ultraviolet divergence of 
the gap integrals. 

Clearly, it is essential to correctly compute the magnitude of the gap in QCD. Studies along this line are 
in progress fl^ , |il| . 
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APPENDIX A: THE MEAN-FIELD APPROXIMATION 



Consider a statistical mechanical system at finite temperature T and chemical potential [i where fermions 
interact via exchange of an TV-component bosonic field. For N = 1 the Lagrangian is identical with eq. ([!]). 
The grand partition function of this system is 

Z =N [ Vi> T>tpT>cf> exp {1$, ip, <()}} , (Ala) 



I[, :,:,,} = I [ ij(x) [G+] \ X ,y)^{y)-l £ <t> a (x) D~^(x, y) 4> h (y) 

a,b=l 

N 



E 

a=l 



^ grP(x)T a i;(x)r(x) . (Alb) 



Here, 

[G±Y\x,y) = -i[i 1 -d±n l0 ~m] 5^\x - y) . (A2) 

D a b{x,y) is the boson propagator, the structure of which need not be further specified at this point. The 
above form includes scalar interactions for N = 1, Ti = 1, vector interactions for N — 4, (ff — g a ~ 1,A1 </> M , 
To = 5a-i,p7 M , D~2 = 9a,-x,ij,[D~ 1 ] >11 ' gb-i,ui v = 0, . . . , 3, and similarly for other interactions. The bosonic 
fields can be formally integrated out with the result 

Z = M' (det D- 1 )' 1/2 J VrfVip exp{/[^,V]} , (A3a) 
I$,il>]= I (ip(x) [G^]~ 1 (x,y)tp(y) + y^^(i)r„i/)(3;) D ab (x,y)ij(y)T b tp(y) J . (A3b) 



x.y 



The last term physically corresponds to the current-current interaction displayed in Fig. ^. Since it is 
biquadratic in the fermion fields, the inte gration over tjj, ip cannot be carried out. In this Ap pendi x, we 
discuss mean- field approximations to (A3a). The idea is to approximate the biquadratic term in ( |A3b ) by a 
bilinear term times a fermion condensate, which then allows for integrating over tp, ip. In principle, one can 
either have an antifermion-fermion condensate or a fermion-fermion condensate, the latter being the case of 
interest in describing the phenomenon of superconductivity. 




FIG. 7. The current- current interaction. 



1. The mean-field approximation for an antifermion-fermion condensate 



In the mean-field approximation one approximates two of the four fermion fields in the last term in (A3b) 



by their expectation value in order to obtain a bilinear form in the fermion fields which allows for integration 



over ip, ip in (A3a). One possibility is to contract ip and ip, as shown in Fig. H (a). The expectation value 



(ip ip) corresponds to an antifermion-fermion condensate. 
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More specifically Ipjj , define 

j a (x) = gtp(x)T a tp(x) , 
and introduce the fluctuation p a (x) of j a (x) around its expectation value (j a (%))■> 

Pa(x) = j a {x) ~ (Ja(x)) . 

Then, assuming D ab (x,y) = D ba (y,x), one has to first order in p a : 

£ E $(x) T a 1>(x) D ab (x, y) ${y) T b ^(y) = ~ v) Mv) 



a , b 



a , 6 



\ E [<?.(*)) ^ Qb (^ y) 0'6(i/)) + 2 ^ W ^"(^ y) (My))} 

a, 6 

J E [-(?-(»)) D ab (x, y) (j b (y)) + 2j a (x) D ab (x, y) (j b (y))} 



a . b 



(A4) 
(A5) 



(A6) 



Inserting this back into (A3b) and integrating over the fermion fields in (A3a) gives the partition function 
in the mean-field approximation for an antifermion-fermion condensate 



Z { p F) = N' (det D 



det 



G 



(FF) 



exp 



Ta iKx)) D ab (x, y) $(y) T b j,(y)) 



x <y o,6 



where the (inverse) fermion propagator is 



G 



(FF) 



(a;,I,) = [G+] '(*,!/) -iff 2 ^*"!/) / a D ab (x , z) (ip(z)r b ip(z)) . 



(A7) 



(A8) 



a , 6 



This last equation is obviously a self-consistency equation (or gap equation, or Schwinger-Dyson equation) 
for the mean-field propagator, since 



(^(z)r b ^(z))^Tr{G+ F> (z,z)r h } 



(A9) 



The explicit evaluation of the partition function is facilitated assuming translational invariance which allows 
to Fourier-transform all quantities to momentum space, for details, see [^lj . 



<¥(y)¥(y)> 




<w (x) \|/(y)> 
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FIG. 8. (a) The mean-field approximation for an antifermion-fermion condensate, obtained by contracting tp(y) 
and ip(y). (b) The mean-field approximation for a fermion-fermion condensate, contracting tp(x) and V>(y). 
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2. The mean-field approximation for a fermion-fermion condensate 



In the previous subsection A 1 , we have discussed the mean-field approximation for an antifcrmion-fcrmion 
condensate. Now we discuss condensation of fermion pairs. Instead of contracting ip and ip, we now contract 
■0 and ip (or equivalently, ip and ip), cf. Fig. || (b). This leads to a fermion-fermion condensate. More 
specifically, introduce the charge-conjugate spinor ipc via (ft), and define f a as 



T a = cv T a c- 1 . 

Then, the four- fermion term in ( A3b| ) is written as 

${x) T a ip(x) ${y) T b ip(y) = - [$c{x) fa ipc{x) tp(y) T b ip(y) + h.c] 

= ~^ Tr [70^(2/, a;) 7o f Q J(x,y)T b + h.c] , 

where we introduced the 4x4 matrix 

Japix, y) = tl)c a (x)$p(y) , Jp a (v> x ) = [ J ai3{ x > y)] f = [7o ip(y)]/3 [$c{x) To] 

Introducing the fluctuation of this matrix around its expectation value, 

p(x,y) = J(x,y) - (J(x,y)) , 



(A10) 



(All) 

(A12) 
(A13) 



and expanding (All) to linear order in p yields 

i>(x) T a ip(x) Tp{y) T b ip(y) 



\ Tr [70 ( J f (y, x)) 70 T a (J (x, y)) T b - 2 70 J f (y, x) 70 T a (J (x, y)) T b + h.> 



(A14) 



The result for the partition function in the mean-field approximation for a fermion-fermion condensate reads 
N' (AetD- l Y 1/2 



7FF) 



with 



where 



exp / ]TTr [ l0 {j\y,x)) l0 t a {J{x,y))Y b +h. C ] D ab (x,y) 

{ Jx -y a.b 

x / VrfVip exp{/[^,-0]} , 



i>{x) [G+] \x, y)^{y) + 2 [i>c{x) & + (x,y)ip(y) + h.< 



x.y 



A+(x, y) = g 2 ^ (Vc(«) $(tf)> T 6 £ afc (x, y) . 

a . b 



(A15) 



(A16) 



(A17) 



This equation uniquely defines the gap equation in the mean-field approximation. In order to solve this 
equation, one has to compute the expectation va lue ( ipc(x) ij)(y)). This is done as follows. The hermitean 
conjugate of the first term in square brackets in ( A16) is 



where 



[ipc(x) A + (x,y)ip(y)] =ip(y) A (y,x)ip c (x) , 



A (y,x) = 70 [A + (x,y)] j 



(A18) 



(A19) 



The first term in ( A16 ) can be rewritten in terms of charge-conjugate spinors as 
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$(x)[G+] 1 (x,y)4>(y)=j c (y)[G Q ] \y,x)^ c {x) 
Introducing the 8-component spinors 



* = 



the action (A16) can be written in compact matrix notation as 

1 



where 



/[*,*] 



s- 1 



<f(x)S- L (x,y)V(y) 



x.y 



A- 
[Go] 



(A20) 



(A21) 



(A22) 



(A23) 



Let us assume translational invariance for the gap matrix, A ± (x,y) = A ± (x — y). Then, the Fourier 
transforms of the fields and \Gq 1 , A^ are 



V 



V 



(A24a) 
(A24b) 
(A24c) 



Note the choice of signs in the exponential factors for ipc and A . From eq. (||) , one would have expected 
the opposite sign. This choice has the consequence 



fe(fc) =C^ T (-k) , ^ c (fc) =^ T (-fc)C , A-(fc)=7°[A+(fc)] t 7 . 
It also ensures that the action is diagonal in momentum space: 



(A25) 



(A26) 



In the conventional approach to superconductivity, the action is only diagonalized after performing a Bogoli- 
ubov transformation. The choice of signs in the Fourier transforms (A24b,A24c) for the charge-conjugate 
spinor and A~ avoids this additional complication. 

In order to complete the calculation of the grand partition function in the mean-field appr oximation, one 
has to perform the Grassmann integration over the fermion fields ip, ip. In the action (A26K however, the 



fields tpc(k ), ipc{k) also enter; these are not independent integration variables on account of ( 3J) . To proceed, 
from ( |A25t ) one derives the identities 

^j(-k)=C^(k) , ^(-fc)=Vo(*)C, (A27) 
and rewrites the integration measure as 

V^Vip = Y[dtp(k) dtp(k) = Y[dtp(k)dip(-k)dtp(k)dilj(-k) = J\fY[dtp(k)dip c (k)dip(k)dtp c (k) , (A28) 



fe>0 



fe>0 



where M is the (irrelevant, since constant) Jacobian from the transformation ( A27 ). Moreover, one can show 
that 



i£*(*)£^*(fc) S £*(*) 



fc>0 



T 



(A29) 
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Then, 



J V^V^ exp {/[§,*]} ee AAdetfex) 



rs- 1 " 




rs- 1 " 




ee Af ^det 


) 


T 


T 





1/2 



The full propagator S(k) is determined from solving 1=5 1 5, with the result 

G+ 



S 



-Gq A+ G 4 



-G+A-G- 



where all functions depend on the 4-momentum k^ and where we have introduced 

G ± = {[G±]~ 1 , S ± eeA^G f A ± . 

The off-diagonal components of S fulfill the identity 



G^ A ± G ± = G T A ± Gfj 1 



(A30) 



(A31) 



(A32) 



(A33) 



This can be prov en directly, or by solving 1 = S S^ 1 . 

From eq. (A31), and from (^ a (k)^'p(k)) = —TS a f3(k) (this identity is proven e.g. in Appendix B of f2l|l 
one obtains 



(Vc(x) i>(y)) = £ £ z- %k - {x ~ v) Go(k) A+(fc) G+(fc) 



V 



Inserting this into (A17) and taking the Fourier transform, one obtains the gap equation 



,T 



A + (fc) =fl 2 -EE f » Dab ( k - «) G o(l) A+ («) G + (q) T b . 



(A34) 



(A35) 



q a, b 



The gap equation has a graphical representation which is derived from Fig. [s| (b) , cf. Fig. ^|. All that changed 
with respect to Fig. || (b) is that the explicit value of (ipc{x) V>(y)) from ( |A31 ) was used. The blob in Fig. 
^ stands for the gap matrix, while the thick (thin) lines represent the full (free) propagator. 

D ab (k-q) 



A + (k) 



A + (q) 



G (q) 

FIG. 9. The gap equation. 



G + (q) 



APPENDIX B: THE STRUCTURE OF A SCALAR GAP MATRIX 



In this Appendix we analyze the Dirac structure of a scalar gap matrix. We find that in general there 
are eight independent gap functions. Fermi statistics imposes a powerful symmetry constraint on these 
functions. With respect to parity, four of the gap functions describe condensation in the channel of even 
parity, the other four in the odd-parity channel. Only fermions with the same helicity condense. With 
respect to chirality, four gap functions describe condensation of pairs with the same, the other four that of 
pairs with the opposite chirality. The chirality and helicity properties of the gap functions suggest a different 
representation scheme for the gap matrix. Within this scheme, we derive that the number of independent 
gap functions reduces to four in the ultrarelativistic limit. 
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1. Dirac structure 



The Fourier transform A(k) of the (translationally invariant) scalar gap matrix A(a; — y) can be expanded 
in the basis of the 16 linearly independent 4x4 matrices. However, for condensation in the J — channel, 
all Lorentz indices have to be contracted with either the 4-momentum, fc M , or, because of the presence of a 
medium, with the respective 4- velocity of this medium, . This leaves the eight matrices 

1, 7-fc, 7-tt, j-kj-u, 75, 757 -k, 757 -u, j 5 j-kj-u. (Bl) 

In the rest frame of the medium, u M = (1, 0), the most general ansatz for A is then 

A = A 1 75 + A 2 7 • k7o75 + A 3 7075 + A 4 + A 5 7 • £70 + A 6 7 • k + A 7 7 • £75 + A 8 7 , (B2) 

where k = k/|k|, and where the gap matrix A as well as the gap functions A„ depend on ko, k. The notation 
follows that of Bailin and Love || . 



2. Symmetry properties 

The antisymmetry of the fermion fields provides a powerful constraint on the gap functions A„. With ([|) 
and C = — C _1 one rewrites 

f ^c(x) A(x - y) i/j(y) = - f ^ T (x) [A(y - x)f $%(y) = f C' 1 [A(y - x)f C^{y), (B3) 

i.e., in Fourier space 

CA(fc)^ 1 = A T (-fc) . (B4) 

Using C 7^ C _1 = — 7J, this implies: 

A n (k) = +Al(-k) , n = l,..., 6 , A n (k) = -A£(-fc) , n = 7,8. (B5) 

If we neglect the energy- momentum dependence of the gap functions A n , as we do in our solution of the 
gap equations, then this equation demands that Ai, . . . , Aq are symmetric matrices in the space of internal 
degrees of freedom, while A7 and Ag are antisymmetric. In QCD, for example, A^ ^ g = A 3 ^ ^, n = 1, . . . , 6, 

while A^ j g — —A J ^ g p n — 7, 8 [for notation see the discussion preceding eq. (Q)]. If the A n do not have 
internal degrees of freedom, then A7 and Ag must vanish, as noted by Bailin and Love ||. We stress, 
however, that A7 and Ag need not vanish when they carry internal degrees of freedom. 



3. Parity 

Under a parity transformation, t — > t, x — * — x, so 3- vectors change their sign. On Dirac spinors, a parity 
transformation is effected by |E2| S(P) = r/p7o, where rjp — ±1 is the intrinsic parity of a particle. Thus 
70 has even parity, S~ 1 (P) 70 S(P) = 70, while 7* has odd parity, S~ 1 (P)j z S(P) 7*- Wave functions 
transform as ip'(t,x) — S(P)ip(t,—x) = rip^ ifj(t,—x), i.e., with (||), ip' c (t,x) — —ipc(t> ~ x ) S(P), which 
shows that the spinor 'ipc has opposite parity from ip. 

We rewrite the action (||) in Fourier space, with the conventions listed in eq. ( A24 ): 

i>]=J2 {^( fc ) [°o ] ~\ k ) V-(fc) + \ [Mk] A(k) rP(k) + h.c] | . (B6) 

Then, the term 

V>c(fc) (A x 75 + A27 • k7o75 + A37075 + A77 • £75) ip(k) (B7) 
represents condensation in the even-parity channel, while the term 

i, c {k) (a 4 + A 5 7 • k 7o + A 6 7 • k + A 8 7o ) ff>(k) (B8) 
represents condensation in the odd-parity channel. 
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4. Helicity 



The spinors ip, %l>, ipci m the effective action (B6) can be decomposed with respect to their helicity. 
The helicity projector is given by 



V ± (k) EE 

and we denote the helicity-projected spinors by 



1 ± 7 5 7o7 • k 



^±(fc)=7>±(k)V(fc) , i>±(k)=rP(k)V±(k) , 
i>c±(k) =P±(k)ipc{k) , ^c± (fc) = $c (k) V± (k) 



(B9) 



(BlOa) 
(BlOb) 



The second equation results from the fact that due to our sign convention in (A24b), in Fourier space 
i)c{k) = V{-k)C. 

The inverse free fermion propagator as well as the gap matrix (B2) commute with the helicity projector, 



[G+] \k),V ± (k)\ =[A(*),P±(k)] = 0. 
This means that the action ( p36| ) does not mix states of different helicity, 



(Bll) 



(B12) 



s=± 



As a consequence, condensation in a scalar ( J = 0) channel can only occur between fermions of the same 
helicity (++ or ). For a physical explanation, see the discussion following eq. (^0|). 



5. Chirality 



The chirality projector is given by 



1±75 



Let us introduce right- and left-handed spinors via 



4>r 



With (|), one then derives 



4>Cr 
IpCr 



IpC V r 



ipce = V r ipc , 



Let us now investigate the effect of the chirality projectors on the gap matrix 
and V"ie = V r ,e, VrVt = one derives 



while 



V r> £ A V r ,t = Ai 75 + A 2 7 • k7o75 + A 4 + A 5 7 • k7 , 



Vi,r A V r j = A 3 7075 + A 6 7 ■ k + A 7 7 • k7 5 + A 8 70 



(B13) 



(B14a) 
(B14b) 



(B15a) 
(B15b) 

3). From 7^75 = —757^ 
(B16) 

(B17) 



This result means that Ai, A2, A4, and A5 are gap functions describing condensation of fermion pairs 
with the same chirality (right-right or left-left), while A3, Ag, A7, and As are gap functions describing 
condensation of fermion pairs with the opposite chirality (right-left and left-right). 

Let us summarize the results of this and the preceding subsections. On ly p articles with the same helicity 
can condense to form a scalar (J = 0) condensate. In the general ansatz (B2), the term 
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Ai 75 + A 2 7 ■ k7o7 5 (B18) 
describes the condensation of fermions with the same chirality in the even-parity channel, while 

A 4 + A 5 7 • k 7o (B19) 
describes condensation of fermions with the same chirality in the odd-parity channel, and 

A 3 7 7 5 + A 7 7-k75 (B20) 
describes condensation of opposite- chirality fermions in the even-parity channel, while 

A 6 7-k + A 87o (B21) 
describes condensation of opposite-chirality fermions in the odd-parity channel. 

6. Quasiprojector representation of the gap matrix 

The properties of the gap matrix with respect to chirality and helicity suggest a somewhat different 



representation. Instead of the eight matrices (Bl), one can use the eight matrices constructed from chirality 
helicity and energy projectors. We use the following projectors onto states of positive and negative energy 
for free particles: 

A ±/ k ^ = ± (7o7 • k + TO70) = 1 ± (/3k7o7 • k + a k 7o) 
2 Efc 2 

where «k = m/E^, /?k = |k|/i?k- These differ by 70 and the normalization of the spinors from the commonly 
used projectors |22|| , but have the advantage that they are regular in the limit m — *■ 0. 
The helicity projectors commute with either the chirality or these energy projectors, 

[V r ,t , V± (k)] = [V± (k) , A± (k)] = , (B23) 

but, for finite m, as massive spinors are not eigenstates of chirality, the energy projectors do not commute 
with the chirality projectors, 

[P r , A± (k)] - T« k 7075 , [Vt , A ± (k)] = ±a k 7o7 5 . (B24) 

Let us introduce the "quasiprojectors" 

Vf t ± (k) = V r ,i V± (k) A± (k) . (B25) 



These eight quantities constitute a basis which is equivalent to (Bl). Therefore, the general gap matrix (B2) 
can alternatively be written as 

A = <p+ + V+ + + (j)+ + V+ + + <i>+_V+_ + 4_V+_ + 07+^7+ + 07+^7f + 07-^7- + . (B26) 

with an obvious notation for the new gap functions 0^£ + - The old gap functions A„ can be expressed in 
terms of the new gap functions 0: 

(B27a) 
(B27b) 
(B27c) 
(B27d) 
(B27e) 
(B27f) 
(B27g) 
(B27h) 



Ax = 


g{(l + /?k) 


w+ - 0/ 


_ + <f, r _ - <f, t+ ] + (1 - 


/3k) [07- - 


07+ + 0r+ - 0£-] } . 


A 2 = 


"{(! + &) 


[-07+ + « 


£- + 07- - 07+] + (! 


- Ac) [07- 


- 07+ - 07+ + 07-] } 


A 3 = 




<t>i- + <t>i 


h + 4>l- + 07+ + 07- - 


- 07+ - 07_ 


] , 


A 4 = 


g {(1 + &) 




_ + 0-_ + 4>j + ] + (1 - 


/3k) [07- + 


07+ + 07+ + 07-] } . 


A 5 = 




* 


+ 07- + ^7+] + C 1 


- /?k) [07- 


+ 07+ - 07+ - 07-] } 


A 6 = 


Y ^ ■+ 


- 07- + <t>i 


f - ^7- + 07+ - 07- - 


- 07 + + 07_ 


] . 


A 7 = 


f [0++-« 


>+_ + <t>i+ 


- 07- - 07f + 07- _ < 


67+ + 07-] 




A 8 = 




£. + 0/+ 


+ 4>i- - 07+ - 07- - ( 


6 7+ - 07-] 
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The new basis ( |B25| ) is complete, 



E££^(k) = i 



(B28) 



h=r,£ s=± e=± 



However, using (B24) we find that the quasiprojectors (B25) are not true projectors: 



v+ ± 


v+ ± = v+ ± + 


-y 7075 VJ ± 


•p+, 1?+ 


= -y 7075^ , 


(B29a) 






y 7o75 


' «± ' fcfc 


, oik 4- 
= K± - y 7o75 "P r ± , 


(B29b) 




v- ± = v r - ± - 


-y 7075 ^± 


^±^± 


= -y 7075 ^± , 


(B29c) 




?Y±7V± = - 


— 7o75 V r ± 


^±^± 


= "P £± + y 7o75 ^r"± > 


(B29d) 




v+ ± v- ± = 


— 7075 V t ± 




= — 2" 7075 Vi± ' 


(B29e) 






«k „_ 

— 7075 V r± 




«k „_ 

= — 2" 7o75 P r ± , 


(B29f) 




v~ ± v+ ± = - 


— 7o75^± 




= -y 7o75 VJ ± , 


(B29g) 




r P7±K± = - 


Y 7075 >V± 




= y 7o75 7V± 


(B29h) 



(the argument k was dropped for convenience). 

For the remainder of this subsection, we consider the ultrarelativistic limit, m = ak = , /3k = 1. In this 
limit, the quasiprojectors become true projectors. Moreover, four of the eight projectors vanish, 



(m = 0) 



(B30) 



This is another way of stating that there are no massless right- (left-) handed fermions with negative (positive) 
helicity and positive energy, and no right- (left-) handed fermions with positive (negative) helicity and 
negative energy. The reduction in the number of projectors can be intuitively understood noting that 
massless fermions can be described in terms of 2-component Weyl spinors which require only the simpler 
algebra of linearly independent 2x2 matrices. 

Another consequence of the ultrarelativistic limit is that one of the three projectors in (B25) is redundant: 



V+ + j_(k) = T r jV±(k) = V Tit A+(k) = V ± (k) A+(k) , 
V-_ /+ (k) = V r jV T (k) = P rA A-(k) ee V T (k) A-(k) . 



(B31a) 
(B31b) 



Thus, we could use any two of the three projectors for chirality, helicity, and energy to construct (B25). We 
keep all three indices, however, because i t fac ilitates the physical interpretation of our results. 

When to = 0, as the four projectors ( B30 ) vanish, the eight independent gap functions reduce to four, 
and eq. (B26) becomes 



v: 



"T— ' r — 



Using eq. ( B27 ) we obtain for the gap functions A n : 







f 


4>t- + 4>7- - <t>7+] 


A 2 


4h 


' } r+ 


+■ 4>~l_ + - (f>j + 


A 4 




f + 


4>t- + 4>7- + <t>7+] 


A 5 


4h 




- 4>~l_ + + cj>j + 


A 3 


= A 6 = 


A 7 


= A 8 = . 



(B32) 

(B33a) 

(B33b) 

(B33c) 

(B33d) 
(B33e) 
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The use of the projectors makes it clear that there are only four condensates for massless fermions. This 
was not apparent previousl y || . 

Another consequence of (B32) is that the action (B6) decomposes into four parts; with V^± = ^P^i^^i 



(B34a) 
(B34b) 



where h = r, £, s = ±, e = ±. We draw the important conclusion that in the ultrarelativistic limit, there is 
condensation only of fermions with the same helicity and the same chirality. 
In the scalar model we find that (j>J, — ~4>t_ and cf)~_ = — 07, ; thus 



A! = \ $ 



A 4 







(B35) 



From eqs. (B18), (BIS) we conclude that condensation occurs only in the + channel, and not the channel. 



APPENDIX C: RENORMALIZING F± 



At T = 0, we rewrite the function Tr: as 



T ± _ 9 

~ T 



d 3 q 

(2tt) 3 



D(k - q, M„) \ - D{k - q, 0) -!- 



with the counter term 



d 3 q , , 1 

^ (k " q ' 0) M 



(Cl) 



(C2) 



We assume that the renormalized boson mass vanishes at zero density and temperature, and so in the counter 
term we take the boso n to have zero mass as well. 

The first integral in (Cl) is now ultraviolet-finite. The counter term (C2) can be computed via dimensional 
regularization. First rewrite (|C2]), using p3j 



i«V> T(a)r((3)J 



dx ■ 



+ b(l-x)] a +P 



and shifting q — > q + krr, into 



CT & Jo * VT= HJ (2^) 3 [q 2 +k2x(l 



a;)] 3 / 2 ' 



Now compute the last integral in d = 3 — e dimensions. Note that this implies g 2 
dimensions of energy. The standard formula [E3| 



dx ar 



1 



(m)r( £ 



g+1 
2 



[x 2 + A/ 2 ]" 2 r(o) [M 2 ]"-^ 1 )/ 2 



leads to 



CT 



16tt 2 



2 , k 2 



(C3) 
(C4) 

g 2 A 3 ~ d , where A has 
(C5) 

(C6) 



where Ae is the renormalization scale, A = Ay'ie 1 ^/ 2 , 7 being the Euler-Mascheroni constant. As usual, 
the 1/e term in (C6) is discarded. 
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